ohibc logo
OHI British Columbia | OHI Science | Citation policy

knitr::opts_chunk$set(fig.width = 6, fig.height = 4, fig.path = 'Figs/',
                      echo = FALSE, message = FALSE, warning = FALSE)

dir_git     <- path.expand('~/github/ohibc')
source(file.path(dir_git, 'src/R/common.R'))  
dir_spatial <- file.path(dir_git, 'prep/spatial')  
dir_anx     <- file.path(dir_M, 'git-annex/bcprep')


### goal specific folders and info
goal      <- 'fis'
scenario  <- 'v2017'
dir_goal  <- file.path(dir_git, 'prep', goal, scenario)
dir_goal_anx <- file.path(dir_anx, goal, scenario)

### provenance tracking
library(provRmd); prov_setup()

### Kobe plot functions
source(file.path(dir_goal, 'kobe_fxns.R'))

### set up proj4string options: BC Albers and WGS84
p4s_wgs84 <- '+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0'
p4s_bcalb <- '+proj=aea +lat_1=50 +lat_2=58.5 +lat_0=45 +lon_0=-126 +x_0=1000000 +y_0=0 +datum=NAD83 +units=m +no_defs +ellps=GRS80 +towgs84=0,0,0'

1 Summary

Process RAM database for British Columbia fisheries stock status and harvest levels. Each RAM stock will be filtered to the RAM stock ID and RAM area (to separate region-specific stocks), and then the stock status parameters B/Bmsy and F/Fmsy for that stock for each year reported.

Total catch data for each stock for each year will be used to determine final weightings. Since the distribution of catch data is not differentiated within the reporting region, we will use SAUP and/or DFO spatial catch data to determine the distribution and allocate the total catch proportionally.


2 Data Source

Reference: RAM

Downloaded: 2016

Description: Fishery status (B/Bmsy and F/Fmsy) for stocks in BC waters

Time range: 1927-2015 (varies by stock)

Format: .RData file


3 Methods

Explore RAM: Preliminary RAM database v. 3.8. Load from .RData file, including time series and metadata. This initial chunk explores the variables included in the data.

3.1 BC-specific stocks

From global RAM data, filter to stocks whose region is ‘Canada West Coast’, and save to the ‘int’ folder in GitHub.

3.1.1 Examine time series of stocks

For each BC-specific stock, determine time series of B/Bmsy, F/Fmsy, and catch.

3.2 Plot stock status against modified Kobe

3.3 Associate stocks with spatial boundaries

Spatial boundaries for RAM stocks are derived from shapefiles provided by DFO.

  • From dfo_khunter/d2017/boundaries.gdb:
    • QCI aka Haida Gwaii, WCVANI, QCS, SoG, CC
  • From dfo_khunter/d2017/shrimp_trawl_areas/DFO_BC_SHRIMP_TRAWL_MGMT_AREA.shp:
    • SMA-14, SMA-FR, SMA-PRD
  • From dfo_khunter/d2016/management_boundaries/groundfish_mgmt_areas/MajorMinor_GroundfishManagementAreas.shp:
    • 5AB
  • BCW and PCOAST are union of all (areas 1-5) from boundaries.gdb
  • WCVANI-SOG-QCI is the union of those areas (1, 4, 5)
  • HS will be the area between Haida Gwaii island and the North Coast, bounded on the south by Area 5AB and bounded on the north at the level of the northeastern-most point of Haida Gwaii

FRDU is Fraser River Designatable Unit - since eulachon are harvested when they return to the mouth of the Fraser River, we will assume that this fishery falls entirely within Strait of Georgia region. If SAUP data are available, perhaps we can use those instead.

## OGR data source with driver: ESRI Shapefile 
## Source: "/home/shares/ohi/git-annex/bcprep/fis/v2017/stock_boundaries", layer: "stock_boundaries_bcalb"
## with 10 features
## It has 2 fields

4 Provenance

  • Run ID: 10 (4260942); run tag: “standard run”
  • Run elapsed time: 55.733 seconds; run memory usage: 21943.2 MB
  • System info:
    • System: Linux, Release: 4.4.0-57-generic. Machine: x86_64. User: ohara.
    • R version: R version 3.4.0 (2017-04-21), Platform: x86_64-pc-linux-gnu, Running under: Ubuntu 14.04.5 LTS.
    • Attached base packages: stats, graphics, grDevices, utils, datasets, methods, base
    • Other attached packages: provRmd_0.1.1, stringr_1.2.0, RColorBrewer_1.1-2, dplyr_0.5.0, purrr_0.2.2, readr_1.1.0, tidyr_0.6.2, tibble_1.3.0, ggplot2_2.2.1, tidyverse_1.1.1
%3 10 kobe_fxns.R 11 DBdata.RData 2 data_prep_fis_ram.Rmd#bc_specific 11->2 used 12 bc_stocks_ram.csv 9 data_prep_fis_ram.Rmd#stocks_b_bmsy_f_fmsy 12->9 used 13 ram_b_bmsy.csv 6 data_prep_fis_ram.Rmd#plot_stock_timeseries 13->6 used 3 data_prep_fis_ram.Rmd#examine_each_stock 13->3 used 4 data_prep_fis_ram.Rmd#get_dfo_boundary_info 13->4 used 14 ram_f_fmsy.csv 14->6 used 14->3 used 15 ram_catch.csv 15->6 used 16 ram_ids_to_names.csv 16->6 used 16->3 used 17 ram_stock_id_spans.png 18 kobe_bocaccio_bc_all.png 19 kobe_pacific_ocean_perch_qci.png 20 kobe_pacific_ocean_perch_wcvani.png 21 kobe_english_sole_hecate_strait.png 22 kobe_canary_rockfish_wcvani-sog-qci.png 23 kobe_rock_sole_qcs.png 24 kobe_rock_sole_hecate_strait.png 25 kobe_pacific_herring_cc.png 26 kobe_pacific_herring_qci.png 27 kobe_pacific_herring_sog.png 28 kobe_pacific_herring_wcvani.png 29 kobe_lingcod_sog.png 30 kobe_pacific_cod_qcs.png 31 kobe_pacific_cod_hecate_strait.png 32 kobe_pacific_cod_wcvani.png 33 kobe_sablefish_bc_all.png 34 kobe_northern_shrimp_sma_fr.png 35 kobe_northern_shrimp_sma_prd.png 36 kobe_sidestripe_shrimp_sma_prd.png 37 kobe_sidestripe_shrimp_sma_fr.png 38 kobe_northern_shrimp_sma_14.png 39 kobe_sidestripe_shrimp_sma_14.png 40 kobe_eulachon_fraser_river_du.png 41 ram_area_ids.csv 41->4 used 42 ram_stock_to_area.csv 43 ohibc_rgn_to_ram_area.csv 1 data_prep_fis_ram.Rmd 8 data_prep_fis_ram.Rmd#setup 1->8 wasExecutedBy 1->2 wasExecutedBy 1->9 wasExecutedBy 1->6 wasExecutedBy 1->3 wasExecutedBy 1->4 wasExecutedBy 5 data_prep_fis_ram.Rmd#match_dfo_areas_to_ohibc_rgns 1->5 wasExecutedBy 7 data_prep_fis_ram.Rmd#provenance 1->7 wasExecutedBy 8->10 wasExecutedBy 2->12 wasGeneratedBy 9->13 wasGeneratedBy 9->14 wasGeneratedBy 9->15 wasGeneratedBy 6->17 wasGeneratedBy 3->18 wasGeneratedBy 3->19 wasGeneratedBy 3->20 wasGeneratedBy 3->21 wasGeneratedBy 3->22 wasGeneratedBy 3->23 wasGeneratedBy 3->24 wasGeneratedBy 3->25 wasGeneratedBy 3->26 wasGeneratedBy 3->27 wasGeneratedBy 3->28 wasGeneratedBy 3->29 wasGeneratedBy 3->30 wasGeneratedBy 3->31 wasGeneratedBy 3->32 wasGeneratedBy 3->33 wasGeneratedBy 3->34 wasGeneratedBy 3->35 wasGeneratedBy 3->36 wasGeneratedBy 3->37 wasGeneratedBy 3->38 wasGeneratedBy 3->39 wasGeneratedBy 3->40 wasGeneratedBy 4->42 wasGeneratedBy 5->43 wasGeneratedBy